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1.  INTRODUCTION 


Numerical  simulation  of  the  ram  accelerator  has  been  an  ongoing  research  project  at 
the  U.S.  Army  Research  Laboratory  (ARL)  since  1990  (Nusca  1990,  i99L  1992).  Computa¬ 
tional  fluid  dynamics  solutions  of  the  Navier-Stokes  equations  have  been  applied  to  turbulent 
non-reacting  and  reacting  in-bore  fiowfields  for  a  ram  accelerator  projectile  launch  system. 
Numerical  simulations  reveal  in-bore  flowfield  details  and  provide  comparison  with  measured 
launch  tube  wall  pressures  and  projectile  thrust  as  a  function  of  velocity.  Numerical  simula¬ 
tions  are  also  used  to  investigate  the  ARL  120mm  (bore  diameter)  ram  accelerator  system. 
These  simulations  are  used  as  a  system  design  aid  and  as  a  means  by  which  geometric  and 
fluid  dynamic  scaling  phenomena  are  investigated. 

For  non-reacting  flow  simulations  the  Rockwell  Science  Center  USA-PG  (Unified  .Solution 
Algorithm  Perfect  Gas)  code  is  used  (Chakravarthy  1985,1989).  In  this  code  lire  Navier- 
Stokes  equations  are  cast  in  conservation  form  and  converted  to  algebraic  equations  using 
upwind  and  central  finite  differences  and  finite- volume  formulations.  The  equations  are 
solved  using  a  second-order  TVD  (total  variation  diminishing)  scheme.  The  code  is  capable 
of  simulating  mixed  subsonic/supersonic  fiowfields. 

For  finite-rate  reacting  flow  simulations  the  ARL-RAMCOMB  (H.AMjet  CO.MBustion) 
code  (Nusca,  May  1990)  has  been  used.  The  steady  2D/axisymmetric  Navier-Stokes  equa¬ 
tions  are  written  in  stream  function-vorticity  form  and  solved  using  a  Gauss-Seidel  relaxation 
scheme.  These  equations  include  conservation  of  chemical  spt^cies  and  reacting  flow  source 
terms  in  the  energy  equation.  Both  one-step/globcd  and  three-step  finite-rate  reaction  mech¬ 
anisms  for  fuel  (CIG)  and  oxidizer  (O2)  are  considered.  In  addition,  the  latest  versions  of 
the  USA-serics  codes  are  being  utilized  (Ota  I9S8,  Paianiswarny  1989).  These  codes  include 
both  equilibrium  (US,A-EC)  and  finite-rate  (USA-RG)  chemistry  which  is  fully  coupled  to 
the  gasdynamics  and  are  capable  of  time-accurate  simulations.  Comparison  of  results  from 
the  USA-RG  and  RAMCOMB  codes  will  be  reported  as  well  as  comparisons  between  the 
USA-EC  code  and  fi  nite-rate  computations. 

Results  for  the  ram  accelerator  projectile  excluding  the  obturator  have  been  published 
(Nusca  1990,1991,1992).  These  numerical  simulations  used  non-reacting  or  reacting  flow 
codes  (in  the  case  of  the  reacting  flow  previous  efforts  used  a  coupling  procedure  for  tlie 
USA-PG  and  ARL-RAMCOMB  codes)  to  simulate  the  effects  of  acceleration  on  the  projectile 
flowfield.  Numerical  simulations  provided  a  possible  explanation  of  the  thrust /velocity  curve 
as  well  as  data  on  the  fluid  dynanne  and  reacting  flow  scaling  effects.  'Phese  results  are  not 
repeated  in  the  present  report  and  the  reader  is  referred  to  Nusca  (Nov.  1991)  for  dotaiis. 
It  is  the  purpose  of  the  present  report  to  investigate  the  projectile/obturator  separation 
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and  combustion  starting  process  using  non-reacting  flow  codes  (USA-PG)  as  well  as  a  fully 
coupled  finite-rate  chemistry  code  (USA-RG).  Whereas  previous  simulations  have  been  for 
steady  and  quasi-steady  flows,  the  present  report  focuses  on  time-accurate  |)rediciions  of 
time-dependent  flow  phenomena. 

The  ram  accelerator  projectile  geometry  u.sed  in  these  studies  represents  a  sim]>lification 
of  the  actual  configuration  used  for  lest  firings.  The  actual  projectile  includes  a  set  of  four 
bore-riding  fins  that  extend  from  the  point  of  maximum  projectile  diameter  to  the  projectile 
base  and  span  the  area  between  the  projectile  and  the  launch  tube.  Exclusion  of  these  fins 
permits  a  2D/axisymmetric  calculation  and  results  in  a  significant  computer  time  savings 
over  the  full  3D  numerical  simulation.  Such  3D  simulations  have  been  performed  with  the 
present  codes. 

2.  BACKGROUND 

The  ram  accelerator  technique  was  first  investigated  via  experimental  test  firings  at  the 
University  of  Washington  (UW)  (Hertzberg  1988,  Bruckner  1988/1991,  Knowlen  1992)  The 
UW  ram  accelerator  facility  uses  a  light  gas  gun  (e.g.  helium  driver  gas)  to  accelerate 
projectiles  up  to  1300  m/s.  The  muzzle  is  connected  to  a  perforated  tube  and  evacuation 
tank  which  serve  as  a  dump  for  the  driver  gas  prior  to  entrance  into  the  16m  long  ram 
accelerator  tube.  The  accelerator  tube,  which  can  be  divided  into  sections  separated  by 
diaphragms,  is  filled  with  a  pressurized  fuel/oxidizer  mixture  and  is  instrumented  at  40  axial 
locations.  Instrumentation  consists  of  pressure  transducers,  fiber-optic  light  guides,  and 
magnetic  transducers.  Thin  magnetic  sheets  are  mounted  in  the  nose-body  joint  and  in  the 
base  of  the  projectile  (see  Figure  1).  When  the  projectile  passes  electromagnetic  transducers 
on  the  accelerator  tube  these  magnets  induce  signals  that  are  used  to  determine  the  distance¬ 
time  history  (i.e.  velocity)  of  the  projectile.  Projectile  thrust  is  derived  from  the  velocity 
history.  The  ARL  120mm  (bore  diameter)  facility  is  described  by  Kruezynski  (1992). 

The  projectile  consists  of  an  axisymrnetric  cone-boattail  body  (blunt  base)  with  stabiliz¬ 
ing  fins  to  center  it  along  the  launch  tube  axis  (see  Figure  1).  The  projectile  is  shaped  like 
the  centerbody  of  a  ramjet  engine  and  is  injected  into  a  stationary  tube  filled  with  a  pres¬ 
surized  gas  mixture  of  hydrocarbon  fuel  (e.g.  CH4),  oxidizer  and  diluents  such  as  CO2,  N2, 
He,  and  Ar.  There  is  no  propellant  on  board  the  projectile.  The  tube  resembles  the  outer 
cowling  of  a  conventional  ramjet  engine.  When  the  injection  velocity  is  greater  than  the 
sound  speed  of  the  gas,  a  strong  oblique  shock  system  develops  on  the  projectile  which  sus¬ 
tains  combustion  around  the  projectile.  In  this  way  the  energy  release  process  travels  along 
with  the  projectile.  Thrust  is  generated  by  the  action  of  high  pressure  reacting  gases  on  the 
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rear  part  of  the  projectile.  Various  combustion  ignition  fneciiani.,;;;s  I,..ve  been  investigated 
including  a  perforated  pusher  sabot.  The  initial  gas  pressure,  fuel/oxidizer  coiupositioie  and 
.sound  speed  can  be  selected  to  achieve  the  desired  acceleration  and  projectile  \elocitv  at  tube 
exit.  Diluents  are  used  to  tailor  tlie  acoustic  speed  of  the  mixture  so  that  the  initial  Mach 
number  of  the  projectile  exceeds  the  minimum  required  (~  2.8)  to  start  the  diffus<'r  (i.e. 
projectilc/tube  clearance  at  maximum  projectile  body  diameter)  atid  tailors  the  h<'at  rehutse 
of  combustion  to  a  level  that  stabilizes  the  shock  system  on  the  projectile  })odv.  Excessive 
diluent  results  in  low  projectile  acceleration  levels,  wher<-a.s  insufheient  diluent  conc<*ntration 
can  cause  pre-ignition  on  the  projectile  forebody  and  deceleration.  In  addition,  fne!  lean 
mixtures  can  result  in  ignition  on  the  projectile  forebody. 

The  total  force  on  the  projectile  is  composed  of  the  drag  force  oti  the  forebodv  and  a 
thrust  produced  by  the  high  pressure  combustion  products  on  the  projectile  afterbodv  and 
base.  Normally  a  net  thrust  is  obtained  since  the  pressure  of  the  combustion  j)roducts  is 
higher  than  that  of  the  compressed  gases  downstream  of  the  nose  shock.  Cas  mixtures  with 
higher  heat  of  reaction  yield  greater  net  thrust.  Significant  combustion  must  occur  only  on 
the  afterbody  of  the  projectile  in  order  to  maximize  thrust  since  combustion  on  the  conical 
nose  contril)ut('s  to  the  drag  force.  As  the  combustion  moves  farther  behind  the  cone-boattail 
junction  on  the  projectile  (throat),  smaller  fractions  of  the  afterbody  surface  area  are  expo.sed 
to  high  pressure  gases  resulting  in  reduced  thrust. 

Several  modes  of  ram  accelerator  propulsion  iiave  been  investigated  cxpcrinientally  (Bruck¬ 
ner  1988).  The  first  mode  api)lies  to  projectile  operating  velocities  below  the  C'hapman- 
Jouguet  (CJ)  detonation  speed  of  the  propellant  mixture  (Hertzberg  1988)  which  typirallv 
ranges  in  Mach  numl>er  from  2.5  to  4  for  hydrocarbon  fuels.  In  this  mode,  the  thrust  is 
provided  by  the  high  pressure  jnojectile  base  pressure  resulting  from  a  norma!  shock  system 
stabilized  on  the  body  by  theorized  thermal  choking  of  the  reacting  flow  at  full  lube  area 
behind  the  projectile  (see  Figure  2).  One  theoretical  model  of  the  thermally  choked  mode 
predicts  that  the  normal  shock  recedes  along  the  body  as  the  projectile  velocity  increases 
(Bruckner  1988).  If  the  projeCvile  afterbody  were  to  taper  to  a  point  and  the  flow  were  invis- 
cid.  this  normal  shock  would  gradually  fail  back  to  the  full  lube  area.  .A  normal  shock  in  a 
constant  area  duct,  followed  by  heat  addition  and  thermal  choking  in  steady  flow,  constitutes 
a  CJ  detonation  wave.  Thus,  theory  predicts  that  the  thrust  goes  to  zero  as  the  projectile 
velocity  approaches  the  CJ  detonation  speed  of  a  particular  propellant  mixture  (Bruckner 
1988).  For  projectile  speeds  as  high  as  85%  of  the  CJ  detonation  speed  of  the  mixture.  F\V 
has  observed  that  the  llirnst  as  a  function  of  Mach  number  is  indeed  accurately  predicted 
by  the  one-dimensional  theoretical  model  of  thermally  choked  mode  (Bruckner  19SS).  At 
projectile  velocities  above  85%)  of  the  CJ  detonation  speed,  however,  the  thrust  typically 
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begins  to  exceed  that  predicted  by  the  tiieoretical  model,  reaching  a  ininiimim  at  velocities 
near  95%  of  the  CJ  detonation  speed  of  the  mixture,  and  then  tends  to  increase  with  ve¬ 
locity.  This  trend  coincides  with  the  experimentally  observed  combustion  occurring  on  the 
projectile  body.  Combustion  on  the  surface  of  the  projectile,  as  opposed  to  solely  in  the 
projectile’s  wake  was  proposed  by  the  present  author  based  on  numerical  simulations  (Nusca 
1991).  Experiments  have  shown  that  in  this  second  mode  called  "transdetonalive''  (‘  :tically 
Mach  number  4  to  G),  the  projectile  can  accelerate  smoothly  at  or  slightly  above  the  C.l  det¬ 
onation  speed  of  the  mixture.  For  operating  at  higher  Mach  numbers,  a  “superdelonative" 
mode  has  been  investigated  (Bruckner  1988).  This  mode  operates  at  velocities  greater  than 
the  CJ  detonation  speed  of  the  mixture.  UW  has  proposed  that  this  supersonic  combustion 
process  involves  shock-induced  combustion,  wherein  the  mixture  is  ignited  by  one  of  sev¬ 
eral  reflected  oblique  shock  w’aves  on  the  body.  The  super.sonic  heat  release  raises  the  gas 
pressure  on  the  afterbody  of  the  projectile,  resulting  in  thrust  as  reacted  propellant  expands 
back  to  full  tube  area. 

Typically,  beyond  the  superdetonative  mode,  measured  projectile  thrust  rajridly  decrcase.s 
perhaps  due  to  extensive  combustion  on  the  forebody  of  the  projectile.  Each  of  the  propulsion 
modes  described  above  and  observed  in  experiments  at  the  UW  has  also  been  observed  in 
numerical  simulations  (Nusca  1990,1991,1992). 


3.  NON-REACTING  FLOW  SIMULATION 


The  Reynolds- Averaged  Navier-Stokes  (RANS)  equations  for  2D/axisymmetric  flow  are 
written  in  the  following  conservation  form.  The  dependent  variables  u.  ?\  and  c  are  mass- 
averaged. 
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where  o  =  1  for  axisymmetric  flow  and  0  for  two-dimensional  flow.  Normal  "s  (a),  shear 
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stress  {Txr)^  heat  transfer  (9),  and  internal  energy  (c)  are  defined  in  Nusca  (Nov.  1991).  The 
flow  medium  is  assumed  to  be  a  perfect  gas  satisfying  the  equation  of  state, 

P  =  P^T  (2) 

The  following  power  law  was  used  to  relate  molecular  viscosity  to  temperature  (Mazor  1985): 


where  po  =  0.1716  MPa,  To  =  491.6  R,  and  n  =  0.64874.  The  laminar  and  turbulent  Prandtl 
numbers,  Pr  and  Prj,  were  assumed  constant  with  values  of  0.72  and  0.9  respectively.  The 
ratio  of  specific  heats,  7,  was  also  assumed  constant.  The  specific  heat  rapacities  at  constant 
volume  and  pressure,  c„  and  Cp,  are  related  as  7  =  Cp/c^.  and  3?  =  (7  -  l)Cp/7. 

Assuming  a  time-invariant  grid  and  using  the  transformation  of  coordinates  implied  by 
T  =  =  Tj{x,y),  Equation  1  can  be  recast  into  conservation  form  where  ^ 

and  j;  are  the  new  independent  variables  and  .t^,  Xp,  t/^.  and  t/p  are  the  four  transformation 
coefficients  obtained  numerically  from  the  mapping  procedure  ( Chakra varthy  1985). 

^  +  Gjy  -  ///y]  =  0  (4) 

The  “Area”  in  Equation  4  denotes  the  area  of  the  finite  volume  cell  under  consideration  at 
the  time  of  discretization  of  the  equations. 

The  shock/boundary-layer  interference  flowfield  between  projectile  and  launch  tube  as 
well  as  the  projectile  wake  can  include  regions  of  recirculating  flow.  Modeling  of  these  regions 
can  be  critical  to  the  overall  flowfield  solution  quality.  However,  most  existing  turbulence 
models  either  do  not  treat  such  regions  or  do  so  in  a  semi-empirical  fashion  that  is  frequently 
inadequate.  To  improve  the  predictive  capability  of  separated  flows  using  R.ANS  codes  a  new 
turbulence  model  has  been  developed  (Goldberg  1986).  The  model  is  based  on  experimental 
observations  of  detached  flows.  The  model  prescribes  turbulence  kinetic  energy  [k)  and  dis¬ 
sipation  (e)  analytically  within  backflows.  A  formula  for  the  eddy  viscosity  (/q)  distribution 
within  backflows  is  derived  and  used  for  the  RANS  equations  when  the  calculations  are  done 
inside  separation  bubbles.  Outside  of  them,  another  turbulence  model  (Baldvvin-Lomax, 
1978)  supplies  the  values  of  eddy  viscosity.  While  the  Baldwin-Lomax  turbulence  model  is 
used  to  detect  flow  separation  and  to  initiate  application  of  the  backflow  model,  the  latter 
model  can  relocate  the  separation  point. 
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4.  REACTING  FLOW  SIMULATIONS 


4.1  Chemical  Reactions.  Systems  of  chemical  reactions  are  usually  complex.  A 
hydrocarbon  fuel,  for  example,  may  contain  hundreds  of  distinct  chemical  species.  During 
combustion,  numerous  short-lived  intermediate  species  are  formed  in  addition  to  the  final 
products  that  include  H2O  and  CO2.  Any  increase  in  computational  accuracy  by  including 
all  intermediate  species  and  reaction  steps  within  the  framework  of  a  calculation  could  be 
defeated  by  uncertainties  in  the  thermodynamic,  transport,  and  chemical-kinetic  properties 
of  these  species.  In  addition,  the  computational  costs  of  a  given  reaction  mechanism  depend 
primarily  on  the  number  of  chemical  species  included,  rather  than  on  the  number  of  reactions 
(Westbrook  1981).  It  is  possible  to  make  useful  predictions  of  combustion  processes  by 
confining  attention  to  a  postulated  global  reaction  between  hydrocarbon  fuel  and  oxidizer 
that  neglects  intermediate  steps  (Westbrook  1981). 

Fuel -f- <502 -f- Diluents  =>  (1 -f- ^)(C02  +  H2O) -f  Diluents  (5) 

where  <5  is  the  stoichiometric  oxidizer/fuel  ratio.  Several  limitations  of  a  global  reaction  model 
should  be  noted  (Rogers  1983).  During  the  ignition  delay  period,  when  the  complete  reaction 
mechanism  will  predict  the  gradual  buildup  of  free  radicals  with  little  or  no  perceptible 
temperature  change,  a  global  model  will  indicate  an  immediate  exothermic  (or  endothermic) 
reaction  with  resulting  temperature  increase  (or  decrease).  Secondly,  the  adiabatic  flame 
temperature  predicted  by  the  global  model  will  be  higher  than  for  the  complete  multi-step 
reaction  mechanism.  This  results  from  the  fact  that  a  global  model  does  not  include  all 
the  product  species  actually  present  in  the  reaction.  If  very  low  pressures  are  avoided,  the 
discrepancy  in  flame  temperature  will  not  be  significant  (Rogers  1983).  This  overestimate  of 
the  adiabatic  flame  temperature  increases  with  larger  equivalence  ratio  (fuel/oxidizer)  and 
is  directly  related  to  the  amounts  of  CO  and  H2  in  the  reaction  products  (Westbrook  1981). 
In  addition,  flame  speed  tends  to  decrease  with  increasing  pressure  for  most  hydrocarbon-air 
mixtures;  thus,  it  may  not  be  possible  to  reproduce  both  high  and  low  pressure  ranges  with 
a  single  reaction  rate  expression  (Westbrook  1981). 

Hydrocarbon  reactions  arc  commonly  used  for  ram  accelerator  testing  at  the  University 
of  Washington  (Bruckner  1988)  for  example, 

|cH4  +  2O2 -h  6N2  (6) 

The  fuel  equivalence  ratio.  is  a  measure  of  the  amount  of  fuel  available  for  combustion. 
When  =  1,  the  proportions  of  fuel  to  oxygen  are  stoichiometric.  For  <t>  values  smaller 
than  unity  fuel  should  be  completely'  consumed.  For  the  reaction  stated  above  4>  =  2.7 
which  is  considered  slightly  outside  the  range  of  well  understood  CH4/O2  chemical  kinetics 
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(Anderson  1991).  This  reaction  is  one  of  a  general  class  of  hydrocarbon  reactions  investigated 
by  Westbrook  and  Dryer  (1984), 

CaHb0c  +  “(a+^  +  c)02  =»  aC0  +  ^H20  (7) 

where  elemental  and  diluent  species  N2,  C2,  and  H2  are  ignored  thus  (a=l,  b=  l.  c=l), 

CH4  +  2O2  =>  CO2  +  2H2O  (8) 

CH4  is  considered  the  fuel,  O2  the  oxidizer,  and  products  are  CO2  and  H2O.  A  three-step 
reaction  mechanism  has  also  been  proposed  by  Westbrook  (1991)  in  which  CH4  is  oxidized 
to  CO  and  H2  with  subsequent  oxidation  of  CO  and  H2  forming  CO2  and  H2O- 

CH4  +  2O2  =>  C0-f2H2+^02  (9) 

CO +  ^02  =>  CO2 

■2H2  -f  O2  2H2O 


4.2  Reaction  Rates  and  Kinetics.  Conservation  of  chemical  species  can  be  written 
for  each  specie  in  Equations  S  and  9.  Because  a  rate  term  is  included  the  governing  equations 
are  mathematically  stiff.  Stiffness  can  be  defined  as  the  ratio  of  the  largest  to  the  smallest 
time  scale.  In  reacting  flows,  timescales  associated  with  chemistry  tend  to  be  much  smaller 
than  time  scales  of  the  fluid  motion,  sometimes  by  orders  of  magnitude.  In  the  present  code 
the  stiffness  problem  is  resolved  by  uncoupling  the  fluid  dynamics  equations  from  the  species 
conservation  (rate)  equations  for  a  fixed  number  of  numerical  iterations  (or  sub-iterations 
between  time  steps),  then  updating  the  rate  terms.  The  reaction  rate  is  defined  (Westbrook 
1981,  1984)  using  the  Law  of  Mass  .Action  and  an  Arrhenius  expression  for  C,  the  specific 
reaction  rate  constant. 


^  -  / 
=  AT^expf^ 


-Eg 


(10) 


where  and  =  4  (for  Eq.  8),  AT^  is  the  collision  frequency,  the  exponential  term  is 
the  Boltzmann  factor,  and  Eg  is  the  activation  energy.  The  reactant  mass  fraction  terms 
fhe  1-step  reaction  mechanism  (Eq.  8)  are  raised  to  non-stoichiometric 
coefficients  in  order  to  simulate  the  rate  limiting  step  (C0-l-|02  =>  CO2)  of  the  S-step 
mechanism  (Eq.  9).  For  reaction  sets  given  by  Equations  8  and  9,  Westbrook  (1977,  1981, 
1984)  specifies  A,  Eg,  a,b.c  and  d  with  q  =  0  and  ^  =  8.314  J/K-mole  (sec  Table  1). 


Westbrook  and  Dryer  (1981)  used  a  simple  procedure  for  obtaining  the  single-step  reac¬ 
tion  rate  expression  (Eq.  10).  Values  for  the  concentration  exponents  a  and  b  were  held  fixed 
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Table  1.  Reaction  Rate  Equation  Data 


Reaction  j 

Sa  (kcal/mole) 

A  (moles/cm^-s 

)  a 

b 

c 

d 

Eq.  8 

30.0 

8.3x10" 

-.3 

1.3 

0.0 

0.0 

Eq.  9,  step  1 

30.0 

8.3x10" 

1.0 

2.0 

0.0 

0.0 

Eq,  9,  step  2 

44.7 

3.0x10*^ 

0.0 

0.5 

1.0 

0.0 

Eq.  9,  step  3 

37.6 

1.0x10" 

0.0 

1.0 

0.0 

2.0 

at  the  values  stated  above.  They  found  that  the  relation  a  +  b  =  I  was  necessary  in  order  to 
properly  reproduce  the  correct  dependence  of  flame  speed  on  pressure  for  gas  mixtures  where 
4>  is  not  unity.  As  expected  the  flame  speed  depends  strongly  on  a  for  fuel-rich  mixtures. 
For  a  =  —.3  the  fuel  acts  as  an  inhibitor  which  matches  observations  for  methane  ignition 
in  shock  tubes.  The  activation  energy,  Ea,  was  also  held  fixed  and  the  pre-exponential  fac¬ 
tor  A  was  varied  until  the  model  correctly  predicted  me^lsured  flame  speed  for  atmospheric 
pressure  and  4»  =  1.  The  resulting  rate  expression  was  then  used  to  predict  flame  speed  for 
other  values  of  pressure  and  4>.  Each  set  of  rate  expression  parameters  was  evaluated  on  the 
basis  of  how  well  it  reproduced  experimental  data  (flame  speed).  The  results  showed  that 
varying  Ea  from  26-50  kcal/mole  (for  Eq.  8)  effected  the  flame  thickness  alone.  Westbrook 
(1981)  suggests  Ea  =  30  and  48.4  kcal/mole  as  the  most  appropriate  values. 


4.3  Governing  Equations.  The  conservation  equations  for  mass  (global)  and  mo¬ 
mentum  arc  the  same  as  those  given  in  Equation  1.  Additionally,  for  reacting  flows,  species 
ma.ss  conservation  is  given  by, 

-Rj  =  0  (11) 

where  Rj  is  the  reaction  rate  per  unit  volume  for  species  j.  Energy  conservation  for  a 
compressible  reacting,  flow  is  expressed  by  the  same  energy  equation  as  used  for  non-reacting 
flow  (see  Eq.  1)  if  expressed  in  terms  of  the  total  enthalpy  (h), 


d  d 

—  [rpinrij  -f  rj^.,.)  -f  —  [rpwrrij  + 


V- 


pVh  A  +  Jfc  -  (ufv  vTe  -f-  WT^} 


=  0 


(12) 


where  J  is  a  flux  term  for  mass  [Jj  =  (/i<,ff/Re)Vmj),  heat  (J/,  =  (/icff/Pr)cp^ VT),  and 
turbulence  kinetic  energy  {Jk  =  (/iefr/Pr)^^')-  The  mass  fraction  and  molar  specific  enthalpy 
for  species  j  are  iiij  and  hj,  respectively.  Radiation  flux  is  neglected  in  Equation  12. 


In  Equation  12  the  shear  stress  (f)  includes  the  Reynolds  stress  w'ith  an  effective  fluid 
viscosity  expressed  as  the  sum  of  the  molecular  and  turbulent  viscosities,  fi^fj  =  p  +  pf 
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A  Wilke’s  mixing  law  is  used  to  compute  the  chemically  reacting  mixture  \i.scosity.  The 
calorically  perfect  gas  cissumption  can  be  made  when  the  temperature  dependence  of  Cp^  for 
the  reactants  and  products  is  not  well  determined.  The  specific  heat  can  also  be  formulated 
using  an  explicit  temperature  dependence  obtained  from  tabulated  data  (Stull  1971). 

CpJ'^  =  Ai  +  A2T  +  +  AiT^  +  A^T^  (13) 


For  N  species  only  —  1  specie  equations  must  be  solved,  since  the  sum  of  the  mass 
fractions  must  equal  unity.  In  effect  the  global  continuity  equation  is  the  ATh  specie  equation 
since  the  summation  of  all  specie  equations  yields  the  continuity  equation.  The  mixture 
equation  of  state  for  a  thermally  perfect  gas  follow's  from  Dalton’s  Law, 


N 


P^pJtTY, 


m. 


M, 


(14) 


where  3?  =  Adj,  Mj  is  the  molecular  weight  of  species  j\  and  3?  is  the  specific  gas 
constant.  Mixture  temperature  (T)  is  obtained  from  the  definition  of  the  stagnation  enthalpy, 


A' 


h  =  + 

j 


_l_1 

Pr 


Sc 


1  1 


Pr 


+ 


Sc 


J. 

PrJ 


m, 


(15) 


with  V  =  {u^  +  and  V  is  the  magnitude  of  the  turbulent  (fluctuating)  velocity. 

The  Schmidt  number  (Sc)  is  assumed  to  be  unity.  The  Prandtl  number  (Pr)  is  assumed  to 
be  nearly  unity  (.9)  which  is  considered  adequate  for  gaseous  flows  even  with  combustion 
(Bradshaw  1981).  Mixture  viscosity  (p)  is  defined  using  Sutherland’s  expression  (Ames 
Research  Staff  1958)  for  T  <  3400R  and  using  Equation  3  for  higher  temperatures. 


A  two-equation  turbulence  model  has  been  suggested  by  Kim  and  Chung  (1989)  for 
reacting  flows.  This  model  describes  the  turbulence  viscosity  (pt)  as  a  function  of  turbulence 
kinetic  energy  {k)  and  dissipation  rate  (c)  as  pj  =  pC^k^ jc.  A  set  of  partial  differential 
equations  is  written  for  k  and  t. 
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where,  Hk  =  n  =  //  +  Ht/KAk  =  =  1.3,  Ci  =  1.44,(72  =  1.92,  C3  =  .09. 

These  equations  are  solved  along  with  the  Navier-Stokes  equations.  Boundary  conditions 
for  k  and  t  are  fc  =  0, e  =  .056fi{dufdy)^fp  for  solid  walls  and  k  =  10“^l^,t  =  k^-^/L  for 
freestream.  Initial  conditions  are  k  —  koo,(.  =  A-*  ^CT”'^/(.3Tx“*Re“‘^  “' ). 

5.  NUMERICAL  ALGORITHMS 

The  spatial  discretization  technique  for  the  equations  of  motion  must  be  reliable  and  ro¬ 
bust  if  it  is  to  successfully  capture  the  complex  physics  of  in-bore  projectile/launch  tube  in¬ 
teracting  flowfields.  The  TVD  formulation  for  the  convective  terms  along  with  a  special  treat¬ 
ment  of  the  dissipative  terms  provides  an  appropriate  simulation.  In  recent  years,  TVD  for¬ 
mulations  have  been  constructed  for  shock-capturing  finite-difference  methods  (Chakravarthy 
1985).  Near  large  gradients  in  the  solution  (extrema)  TVD  schemes  automatically  reduce  to 
first-order  accurate  discretizations  locally  while  away  from  extrema  they  can  be  constructed 
to  be  of  higher-order  accuracy.  This  local  effect  restricts  the  maximum  global  accuracy  possi¬ 
ble  for  TVD  schemes  to  third  order  for  steady-state  solutions.  TVD  methods  manifest  many 
properties  desirable  in  numerical  solution  procedures.  By  design  they  avoid  numerical  oscil¬ 
lations  and  “expansion  shocks”  while  at  the  same  time  being  of  higher-order  accuracy.  TVD 
formulations  are  also  based  on  the  principle  of  discrete  or  numerical  conservation  which  is 
the  numerical  analog  of  physical  conservation  of  mass,  momentum,  and  energ}*.  Thus,  TVD 
schemes  can  “capture”  flowfield  discontinuities  (e.g.  shock  waves)  with  high  resolution.  At 
a  fundamental  level  they  are  based  on  upwind  schemes;  therefore,  they  closely  simulate  the 
signal  propagation  properties  of  hyperbolic  equations.  Schemes  ba.sed  on  the  TV'D  formula¬ 
tion  are  completely  defined.  In  contrast,  schemes  based  solely  on  central  differences  inv'olve 
global  dissipation  terms  for  stability  and  have  one  or  more  coefficients  that  must  be  judi¬ 
ciously  chosen  to  achieve  desirable  results.  Any  conventional  time  discretization  method 
suitable  for  the  Navier-Stokes  equations  can  be  used  together  with  this  space  discretization 
methodology;  for  example,  approximate  factorization  and  relaxation  techniques. 


5.1  Flowfield  Grids.  Computing  in-bore  projectile  flowfields  is  complicated  by  the 
multi-wall  geometry.  The  ram  accelerator  projectile  consists  of  several  sharp  corners  that 
would  severely  hamper  conventional  grid  generation  schemes  that  require  one  set  of  grid 
lines  to  be  tangent  to  surfaces  and  another  set  to  be  normal  to  them.  This  geometry  is 
more  easily  gridded  by  the  zonal  approach.  The  internal  geometry  of  the  ram  accelerator 
launch  system  in  broken  up  into  three  zones  of  simple  geometric  shape  (zone  1  between  the 
projectile  and  the  launch  tube,  zone  2  in  the  projectile  wake  and  zone  3  aft  of  the  obturator). 
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An  algebraic  grid  is  generated  in  each  zone  with  clustering  near  surfaces  and  other  regions 
in  the  flowfield  where  high  gradients  are  expected.  The  zonal  boundaries  are  transparent  to 
flowfield  phenomena  (e.g.  shock  weaves).  The  actual  grid  used  for  these  computations  (~ 
56,000  nodes)  consisted  of  the  following  dimensions  for  zones  1  and  2/3:  261x41,  741x61. 
The  computational  domain  starts  a  small  distance  forward  of  the  body  and  extends  4.5 
projectile  body  lengths  downstream.  The  computational  results  were  found  to  be  essentially 
independent  of  further  grid  refinement  when  all  other  factors  were  the  same.  A  dense  grid  is 
preferred  for  resolving  flow  details  within  boundary  layers,  to  prevent  shock  smearing,  and 
to  resolve  shock/boundary-layer  interactions. 


6.  RAM  ACCELERATOR  CONFIGURATIONS 

For  the  simulation  of  projectile/obturator  separation,  the  University  of  Washington 
3Smm  system  shot  CSlO  {N2  gas  fill)  wcis  used  (Knowlen  1992).  The  projectile  consisted  of 
a  12.5-degree,  65mm  long  conical  forebody  and  6.25-degree,  50.8mm  long  conical  afterbody. 
The  projectile  base  diameter  was  17.8mm.  The  obturator  was  11.2inm  thick  (12  grams) 
and  consisted  of  5mm  diameter  holes,  19  in  number  or  approximately  33%  porous.  For  the 
reacting  flow  .simulations,  the  projectile  consisted  of  a  10-degree,  83mm  conical  forebody  and 
8-degree,  71mm  conical  afterbody.  The  obturator  was  16mm  thick  and  consisted  of  5mm 
diameter  holes,  19  in  number  or  approximately  33%  porous.  For  comparison  to  the  ARL 
I20min  ram  accelerator  system  these  38mm  system  configurations  were  not  simply  scaled 
but  tiie  actual  120mm  system  geometry  was  used.  The  projectile  consisted  of  a  10-degree, 
261mm  conical  forebody  and  4-degree,  261mm  conical  afterbody.  The  obturator  was  112mm 
thick  (1.1915  kg)  and  was  assumed  to  be  of  the  same  porosity  eis  the  38mm  system.  In  all 
configurations  the  projectile  fins  were  ignored  as  well  as  the  non-porous  obturator  backplate. 

7.  RESULTS 

Simulation  of  the  projectile/obturator  separation  (ignoring  the  separation  of  the  solid 
obturator  backplate)  was  accomplished  by  assuming  that  the  projectile  and  obturator  were 
mated  as  they  entered  the  accelerator  tube  from  the  launch  tube.  Simulation  of  the  launch 
tube  section  of  the  system  is  addressed  by  Nusca  (Oct.  1992).  The  velocity  of  the  projectile 
as  it  traversed  the  accelerator  tube  was  measured  for  the  University  of  Washington’s  shot 
CSIO  and  was  specified  in  the  simulation.  Separation  of  the  obturator  due  to  the  aero¬ 
dynamic  forces  acting  on  both  sides  of  the  disk  was  simulated  by  computing  a  obturator 
force  coefficient  as  the  time-accurate  simulation  was  run.  This  coefficient  along  with  the 
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obturator  mass  was  used  to  update  the  obturator  position  (and  re-grid  zones  2  and  3  of  the 
computational  domain)  with  respect  to  the  projectile,  during  the  run.  The  simulation  was 
stopped  periodically  (corresponding  to  the  location  of  measurement  stations  on  the  Univer¬ 
sity  of  Washington’s  system)  to  examine  the  accelerator  tube  wall  pressure  and  compare  to 
measured  data.  The  results  included  in  this  report  focus  on  the  first  two  stations  (i.e.  up 
to  .751m  from  the  entrance)  although  the  simulation  was  run  for  3.0m  along  the  accelera¬ 
tor  tube  (6m  in  length).  For  CSIO  the  tube  was  filled  with  nitrogen  gas  at  35  atm.  The 
entrance  velocity  of  the  projectile/obturator  was  taken  as  1335  m/s  (from  measurements). 
The  120mm  system  used  for  comparison  has  a  fill  pressure  of  68  atm  and  entrance  velocity  of 
about  1200  m/s.  The  numerical  simulation  was  performed  from  the  fixed  projectile  reference 
frame  with  the  accelerator  tube  moving  at  a  fixed  (or  variable)  velocity  along  with  a  slug 
of  gas  upstream  of  the  projectile.  As  a  result,  the  boundary  layer  on  the  accelerator  tube 
wall  does  not  form  until  the  first  shock  reflection.  The  holes  in  the  obturator  were  treated 
by  assuming  uniform  porosity  of  .33%. 

Figure  3  shows  the  projectile  velocity  versus  distance  along  the  accelerator  tube  as  mea¬ 
sured  at  the  University  of  Washington  (shot  CSlO)  used  for  the  simulations.  The  computed 
obturator  velocity,  initially  the  same  as  the  projectile’s  rapidly  decreased  as  measured  but  at 
values  about  15%  higher.  For  the  120mm  system,  the  obturator  velocity  is  about  10%  lower 
indicating  that  it  separates  more  slowly.  This  is  not  only  due  to  the  maiss  of  the  120mm  sys¬ 
tem  obturator  but  the  lower  pressures  on  the  obturator  surface  (relative  to  those  in  the  38mm 
system).  Scaling  of  the  fluid  dynamics  has  therefore  not  been  completely  accomplished. 

Figure  4  shows  the  computed  obturator  force  coefficient  as  a  function  of  separation  from 
the  projectile.  Initially  the  force  is  high,  but  as  the  obturator  separates  from  the  projectile 
this  force  decreases  (relief  effect)  and  then  increases  with  separation  distance.  The  shielding 
of  the  obturator  from  the  oncoming  flow  (relative  to  the  projectile-based  coordinate  system) 
is  gradually  removed  and  the  drag  rises.  The  force  on  the  obturator  for  the  r20mm  system, 
for  the  same  relative  separation  distance,  is  slightly  smaller  initially  and  slightly  greater  at 
larger  separation.  This  indicates  a  lower  initial  obturator  surface  pressure,  as  mentioned 
previously. 

Figures  5  illustrates  the  flowfields  over  the  projectile/obturator  configurations  at  mea¬ 
surement  stations  1  and  2  (as  well  as  pre-station  1)  in  the  38mm  system.  These  pressure 
contours  highlight  areas  of  large  pressure  gradients  and  show  that  the  projectile  base  flow  is 
effected  by  the  presence  of  the  obturator. 

Figures  6  and  7  show  the  accelerator  tube  wall  pressure  as  measured  in  the  University  of 
Washington’s  .38mm  system  and  computed  in  the  ARL  simulation.  Figure  6  shows  station  1 
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(just  after  entrance  to  the  accelerator  tube)  pressures  where  the  obturator  is  computed  to  be 
located  approximately  78mm  (67%  projectile  length)  behind  the  projectile  base  (note  that 
the  projectile  is  .116m  long).  Pressure  rises  over  the  projectile  are  indicative  of  nose-shock 
reflections  from  the  tube  wall.  These  pressures  are  smaller  than  measured  due  to  the  absence 
of  the  projectile  fins  in  the  simulation.  The  double  pressure  peaks  in  front  of  the  obturator 
location  are  seen  in  the  measured  and  computed  data.  As  indicated  previously  the  obturator 
is  slightly  aft  of  the  measured  location  (higher  velocity,  see  Figure  3)  in  the  simulation.  The 
smaller  pressure  computed  behind  the  obturator  location  is  due  to  the  downstream  boundary 
condition  imposed  (i.e.  tube  fill  pressure  cis  opposed  to  gas  pressure  with  venting  created  by 
the  launch/ vent  tube  which  was  not  part  of  the  present  simulation).  Figure  7  shows  results 
for  station  2  where  the  obturator  is  computed  to  be  approximately  240mm  (200%  projectile 
length).  The  overall  pressure  levels  are  lower  and  the  agreement  between  measured  and 
computed  data  is  better.  The  measured  pressure  rise  at  about  .25m  in  Figure  7  occurs  at 
about  .28m  in  the  simulation.  The  back  pressure  mismatch  is  again  due  to  the  downstream 
pressure  boundary  condition. 

Figures  8  and  9  compare  the  computed  accelerator  tube  wall  pressures  for  the  38mm 
and  120mm  systems.  Due  to  the  different  projectile  lengths,  the  x-axis  has  been  made 
nondimensional.  The  gas  pressure  in  the  3Smm  system  was  35  atm  while  the  120mm  system 
was  charged  to  68  atm.  Since  the  120mm  system  projectile  has  a  different  forebody  angle 
(see  section  “Ram  .Accelerator  Configurations”)  the  pressure  spikes  are  larger  and  at  slightly 
different  locations  for  this  system.  The  pressures  in  the  vicinity  of  the  obturator  are  slightly 
lower  for  the  120mm  system  at  station  1  (Figure  8)  but  slightly  higher  at  station  2  (Figure  9), 
which  is  consistent  with  the  data  in  Figure  4.  As  the  obturator  separates  more  significantly 
from  the  projectile  (Figure  9)  trends  in  pressure  levels  (downstream  of  the  projectile)  are 
comparable  for  the  two  systems.  This  reflects  the  fact  that  some  fluid  dynamic  analysis  has 
been  used  in  the  design  of  the  r20mm  system  projectile  in  order  to  achieve  a  good  scaling 
match  to  the  38mm  system. 

Figure  10  shows  the  results  from  a  projectile/obturator  separation  simulation  for  reacting 
flow,  methane/oxygen/nitrogen  mixture  at  12  atm  fill  pressure  in  the  3Smm  system.  The 
projectile/obturator  geometry  is  slightly  different  (as  detailed  in  the  last  section).  In  this  fig¬ 
ure  the  obturator  is  computed  to  be  approximately  56mm  (37%  projectile  length  of  154mm) 
behind  the  projectile  base.  The  computation  is  for  finite- rate  chemical  kinetics  (three-step 
global  model).  Computed  pressures  before  the  reaction  zone  (before  about  .r2m  in  the  fig¬ 
ure)  compare  well  with  measured  data.  In  the  reaction  zone  pressures  are  slightly  smaller 
(due  to  the  simplified  chemical  kinetics  model)  for  the  computation.  Pressures  near  the 
obturator  location  are  similar  to  measurements  with  smaller  back  pressure  since  the  launch 
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tube  section  has  been  ignored  (see  previous  discussion).  Figure  11  shows  the  projectile  thrust 
in  the  38mm  system  as  a  function  of  velocity.  Overprediction  of  the  thrust  is  expected  to  be 
corrected  when  more  detailed  chemistry  models  (i.e.  more  reaction  steps)  are  included. 

8.  CONCLUSIONS  AND  FUTURE  WORK 

Computational  fluid  dynamics  solutions  of  the  Navier-Stokes  equations  have  been  applied 
to  both  non-reacting  and  reacting  in-bore  flowfields  for  a  ram  accelerator  projectile  launch 
system.  Good  comparison  between  computed  and  measured  pressures  for  non-reacting  flow' 
in  a  38mm  system  was  achieved,  including  the  investigation  of  geometric  scaling  to  a  120mm 
system.  Previous  publications  by  the  ARL  have  focused  on  the  simulation  of  proje?ctile  flow 
phenomena  in  the  accelerator.  The  present  work  represents  numerical  simulation  of  pro¬ 
jectile/obturator  separation  and  the  starting  mechanism  after  projectile/obturator  entrance 
into  the  accelerator  tube.  Good  agreement  with  measured  pressures  during  obturator  sep¬ 
aration  has  been  achieved.  Application  of  more  detailed  reaction  kinetic  models,  as  well  as 
the  investigation  of  the  launch/vent  tube  (pre-accelerator)  sequence,  is  warranted. 
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Figure  3.  Projectile  and  obturator  velocity  along  the  accelerator  tube,  nitrogen  gets  fill 
(35  atm),  38mm  and  120mm  systems. 
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16 


Figure  5.  Computed  pressure  contours  for  nitrogen  gas  fill  (35  atm).  38mm  system, 

three  stations  in  accelerator  tube.  Pre-first  station  with  obturator  at  approx. 
0.14  meters  (computed),  first  station  with  obturator  at  approx.  0.2  meters 
(computed ).second  station  with  obturator  at  approx.  0.36  meters  (computed). 
Axes  are  lengths  in  meters.  Note  magnified  y-scale.  i.e.  angles  on 
projectile  are  not  true. 
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Figure  8.  Computed  tube  wall  pressure  for  nitrogen  gas  fill,  38mm  (35  atm)  and 
120mm  (68  atm)  systems,  first  station  in  accelerator  tube,  obturator  at 
approx.  1.7  L  (computed). 


Figure  9.  Computed  tube  wall  pressure  for  nitrogen  gas  fill,  38mm  (35  atm)  and 

120mm  (68  atm)systems.  second  station  in  accelerator  tube,  obturator  at 
approx.  3.1  L  (computed). 
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Figure  10.  Tube  wall  pressure  for  reacting  2.7CH4-f 202+  5.6N2  gas  fill  (12  atm),  38mm 

systems,  first  station  in  accelerator  tube,  obturator  at  approx.  0.22  m  (measured). 
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Figure  11.  Projectile  thrust  coefficient  vs.  projectile  velocity  for  reacting 

2.7CH4+202+5.6N2  gas  fill  (20  atm),  38mm  system  (A  =  .OOllm  ^). 
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cross  sectional  area  of  launch  tube 
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specific  heat  capacity,  constant  volume 

specific  reaction  rate  constant 

specific  total  internal  energy 

activation  energy 

flux  vectors  (Eq.  1) 

molar  specific  enthalpy 

total  enthalpy 

flow  rate  or  flux 

turbulence  kinetic  energy 

total  body  length 

species  mass  fraction 

molecular  weight 

stoichiometric  coefficient 

number  of  species 

static  pressure 
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heat  transfer  rate 

radial  direction 

reaction  rate  per  unit  volume 
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universal  gas  constant, 
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u  axial  velocity 

V  radial  velocity 

V  magnitude  of  the  local  velocity  vector 

W  dependent  variable  vector  (Eq.  1) 

x,y  cartesian  coordinates 

Greek  Symbols 

a  parameter  (Eq.  1,  10) 

7  ratio  of  specific  heats,  Cp/c^ 

r  diffusion  coefficient 

8  stoichiometric  oxidizer/fuel  mass  ratio 

£  turbulence  dissipation  rate 

Tf  transformed  coordinate 

heat  transfer  coefficient 
/I  molecular  viscosity 

^  transformed  coordinate 

p  density 

a  normal  stress 

T  transformed  time 

Tj-r  shear  stress  tensor 

f  shear  stress  vector 

fuel  equivalence  ratio 

Superscripts 

unit  vector 

total  or  stagnation 

rate 
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Subscripts 


eff  effective 

h  enthalpy 

j  j-th  mixture  component  or  species 

k  turbulence  kinetic  energy 

p  constant  pressure 

r  radial  component  or  radial  direction 

i  turbulence  quantity 

V  constant  volume 

X  i-direction 

z  axial  component 

t  turbulence  dissipation  rate 

T}  7/-direction  transform  coefficient 

0  azimuthal  component 

^  ^-direction  transform  coefficient 

oc  freestream  quantity 
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